Decoding the Mechanism for the Origin of Dark Matter 
in the Early Universe Using LHC Data 
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It is shown that LHC data can allow one to decode the mechanism by which dark matter is gen- 
erated in the early universe in supersymmetric theories. We focus on two of the major mechanisms 
for such generation of dark matter which are known to be the Stau Coannihilation (Stau-Co) where 
the neutralino is typically Bino like and annihilation on the Hyperbolic Branch (HB) where the 
neutralino has a significant Higgsino component. An investigation of how one may discriminate 
between the Stau-Co region and the HB region using LHC data is given for the mSUGRA model. 
The analysis utilizes several signatures including multi leptons, hadronic jets, b-tagging, and miss- 
ing transverse momentum. A study of the SUSY signatures reveals several correlated smoking gun 
signals allowing a clear discrimination between the Stau-Co and the HB regions where dark matter 
in the early universe can originate. 



Introduction. In the near future, data from the 
LHC will be available allowing one to test models of 
physics beyond the Standard Model (SM). Supersymme- 
try (SUSY) and more specifically supergravity grand uni- 
fied models [l|, 0, [1] provide a well motivated framework 
for the exploration of new physics. Thus supergravity 
grand unified models naturally lead to the lightest neu- 
tralino as the lightest SUSY particle, or the LSP, over a 
significant part of the parameter space and with R par- 
ity it is then a candidate for dark matter. An analysis of 
the relic density of the LSP reveals three broad regions 
where the WMAP 0] constraints are satisfied: these in- 
clude (1) the Hyperbolic Branch (HB)[5|, Q where multi 
TeV scalars can appear consistent with small fine tuning 
(this region is alternately referred sometimes as the Focus 
Point region (FP) or as HB/FP), (2) the coannihilation 
regions 0, H, BUS EL (3) the Higgs pole regio nfill (f or 
recent work on these regions(l-3) see [HI, Q Il5l. flgL 
Of these, the stau coannihilation region and the HB re- 
gion are more generic while the pole region (light Higgs 
and the CP odd Higgs A) is more fine tuned. In addi- 
tion there is also the parameter space in the bulk region 
where the relic density is satisfied due to a combination 
of effects. An interesting issue relates to the following: 
to what extent the LHC data will allow one to decode 
the mechanism by which dark matter is generated in the 
early universe. Specifically we will focus on dark matter 
originating in the Stau-Co region or in the HB region to 
answer this question. 

For concreteness we work within the framework of the 
mSUGRA modeld, which is characterized by the pa- 
rameters ttiq , ui\j2 , Aq 7 tan (3 and sign(/i), where mo is 
the universal scalar mass, TO1/2 is the universal gaugino 



mass, Aq is the universal trilinear coupling (all at the 
grand unification scale), and tan/3 = (H2) / (Hi) where 
H2 gives mass to the up quarks and Hi gives mass to 
the down quarks and the leptons, while /1 is the Higgs 
mixing parameter. The parameter space of mSUGRA 
investigated in this work is as in 14| . In this framework 



we will show that the LHC signatures carry sufficient 
information for the discrimination between the Stau-Co 
and the HB regions. 

In the analysis, the sparticle masses and mixings are 
derived from the GUT scale with the SuSpect code [17 1 
coupled to micrOMEGAs [l8[ . We merge the models via 
the SUSY Les Houches Accord format [l!| into Pythia 
[20l | for the computation of SUSY production at the LHC, 
in concert with PGS4 [2l[, to simulate LHC detector ef- 
fects, and obtain the final event record. The models are 
constrained by their ability to properly break electro- 
weak symmety, by sparticle mass limits from LEP and 
Tevatron analyses, FCNC constraints including b — > S7 
and B s — > by the supersymmetric contribution to 

the muon anomalous magnetic moment, and the dou- 
ble sided bound on the relic density (see [H, H2 for 
details). Our post trigger level cuts are as given in 
fl(| and are standard and we list them below: (1) In 
an event, we only select photons, electrons, and muons 
that have transverse momentum Pt(p) > 10 GeV and 
\r](p)\ < 2.4, p = (7, e,fi); (2) For hadronically decay- 
ing tau (jets): P T (r) > 10 GeV and \rj(r)\ < 2.0 are 
selected; (3) For other hadronic jets only those satisfying 
P T (jet) > 60 GeV and \r)(jet)\ < 3.0 are selected; (4) We 
require a large amount of missing transverse momentum, 
P™ ss > 200 GeV; (5) There are at least two jets that 
satisfy the Pt and 77 cuts. Variations on (4) and (5) will 
be discussed later. 
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The SM background was determined through the gen- 
eration of events from QCD multi-jets from light and 
heavy quark flavors, Drell-Yan, single Z/W production 
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along with quarks and gluons (Z/W + jets), and ZZ, 
WZ, WW pair production which give multi-lepton back- 
grounds. We have cross checked our SM backgrounds and 
several other elements of our analysis with simulations 
done by CMS 0, [24[ and the results of these checks are 
in good agreement. We have also found good agreement 
with the SUSY signal analyses of [H, [2J] and the total 
background analysis of [15] under similar cuts. In PGS4 
jets are defined through a cluster-based algorithm which 
has a heavy flavor tagging efficiency based on the param- 
eterizations of the CDF Run 2 tight/loose SECVTX tag- 
ger [25J and is a displaced (secondary) vertex b-tagging 
algorithm which allows detection of b quarks. The b- 
tagging efficiency enters as a product of two polynomials 
each a separate function of \r](jet)\ and Pr(jet). The 
efficiency is maximized in the region |ry(j'et)| < 1 with 
maximal efficiency ef, = (0.4, 0.5) for tight and loose 
tags respectively, and falls off sharply for \r](jet)\ > 1 
with virtually zero efficiency out near |?7(jei)| = 2 and 
Pr(jet) ~ 160 GeV. The analysis of tau decays is done 
using Tauola[26|(For further details regarding PGS4 see 

There are several recent works which discuss dark 
matterSJ^ [3l G|| &l JM J3EJ3I and collider 

physicsHPHIlTMlIfl m, EI El EI Ea ES S3 ■ The 

analysis presented here is very focused in that it con- 
nects LHC signatures directly to the possible origin of 
dark matter, and here we consider two dominant mech- 
anisms, i.e., stau coannihilation and annihilation on the 
HB. We do not go into the details of the relic density 
analysis which are standard. Rather we go directly to 
a discussion of how the experimental data, in particular 
the LHC data, can allow one to decode the mechanism 
for the generation of dark matter. 

Decoding the mechanism for the generation of dark 
matter with LHC data: It is well known that in the stau 
coannihilation region, the neutralino is typically Bino 
like while in the HB region there is a significant Hig- 
gsino component. In the analysis of [Llj it was shown 
that the Stau-Co region is constituted of a collection 
of mass hierarchical patterns, where mSP5 is the domi- 
nant pattern (defined by the mass hierarchy: Xi < "Tl < 
Ir < v T ). Similarly the analysis of [ljj shows that the 
HB region is dominated by the chargino patterns, where 
the chargino xf is the NLSP, and were classified in [l4[ 
as mSPl-mSP4 (where mSPl is the dominant pattern 
defined by the mass hierarchy: X? < xf < X2 < xT)- 
The largest cross sections in the direct detection experi- 
ments arise from a Chargino Wall (CW) [3l| constituted 
uniquely of mSPl. Along the Wall the thermal annihila- 
tion cross sections in the early universe would have arisen 
mostly from X1X1 annihilations into WW, ZZ, ti, bb for 
the Higgsino like LSP. Additionally there are a signif- 
icant number of cases where the annihilation is domi- 



nated by the processes X1X1 — * bb [~ (85 — 90)%] and 
X1X1 t + t~ [~ (5 — 10)%], and these cases typically 
occur for larger neutralino masses which are more Bino 
like, but can also occur for low values of the neutralino 
masses as well on the CW. 

The Chargino Wall referred to above enters impor- 
tantly in the analysis of LHC signatures. It refers to the 
region of HB where the NLSP is the lightest chargino, 
the LSP is mostly higgsino like, and the spin indepen- 
dent cross section is essentially constant O(10 -8 ) pb as 
function of the neutralino mass for neutralino mass in the 
range ~ (80-650) GeV. To explain this feature we begin 
by exhibiting the LSP Xi state in terms of the gaugino- 
Higgsino states so that X? = X = n iB + ri2W 3 + n 3 Hi + 
UaH®- The LSP gaugino-Higgsino content enters impor- 
tantly in the thermal annihilation cross sections that de- 
termine the proper relic density of decoupled neutralinos 
from the epoch of freeze out. It also enters prominently 
in the strength of the scalar neutralino-proton cross sec- 
tion which we now discuss. Thus on the Wall the heavier 
of the CP Higgs mass has a lower limit near 300 GeV and 
more typically it extends into the range of a TeV to sev- 
eral TeV. On the CW one typically has m 2 H ^> m\, and 
sin a « a where a is the Higgs mixing parameter which 
enters in the diagonalization of the CP even Higgs mass 2 
matrix. Further, the sfcrmion poles can be neglected as 
they make a small contribution in this region. Under 
the above limits one finds that the product a x tan (3 is 
essentially a constant, i.e., a x tan/3 ~ — 1. Under this 
circumstance the spin independent cross section in the 
absence of CP phases is given by [H, H^, [5(| 
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x(n 4 + an 3 ) 2 (9/ p + 2/ pG ) 2 



(1) 



Here \x xp is the reduced mass, and f p and f p c are matrix 
elements defined by f p = Ei=«,d, s /f> fpG = (1 - fp) 

where m p f^ = (p\m qi qiqi\p) . The typical ranges for 
rii on the wall are: ni € (.85, .99), n-i <C ri\, and € 
(.1, .6) ~ — 0(114). Using numerical values of / p ,/ p g[48, 
El, H3] one gets erg, (WALL) - 2 x 10~ 8 [pb]. In our 
analysis presented later (see Fig. (5)), however, we have 
implemented the full cross section calculation without 
any of the above approximations. This analysis leads 
to a^ p (WALL) lying in the range ~ (1.5 - 5) x 10~ 8 pb 
while the most recent limits give cr^p ~ 5 x 10- 8 pb[5l| for 
m x « 60 GeV. Thus this region of the parameter space 
is within reach of the current and the next generation of 
dark matter experiments. As noted already, the CW is 
also a very interesting region for LHC signatures. In the 
following we discuss several signatures [listed as (i)-(v)] 
which allow one to discriminate between the Stau-Co and 
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FIG. 1: Left panel: -/Vsusy vs. P™ lss for points in the Stau-Co and HB regions along with the SM background under 
the standard post-trigger level cuts. The Stau-Co and HB model points given here have (mo, miw, Ao, tan/?, sign(u)) as 
(71.5, 348, 334, 10, +) and (1694, 216, -740, 50, +) respectively (all masses are in GeV). The top mass is taken to be 170.9 GeV. 
Right panel: JVsusy vs. (P™ lss ) for each parameter point in the Stau-Co and HB. (P™ lss ) acts as an indicator of Stau-Co and 
HB regions. 



the HB regions over a significant part of the parameter 
space. 

(i) Pj? ISS Distributions: A powerful signature for the 
discrimination of Stau-Co and HB is the total number of 
SUSY events -/Vsusy as a function of the missing trans- 
verse momentum. An analysis of this signature is given 
in Fig.fl]) (left panel) where one finds different geometries 
in these distributions. Here we emphasize that the signif- 
icantly fatter Pj? lss distributions for points in the Stau-Co 
region contrast sharply with the P™ lss distributions from 
points in the HB region which are much thinner [16| as 
exhibited for the case given in Fig. (1) (left panel). Such a 
signature has the interesting feature in that the discovery 
potential is increased over a larger region in the SUGRA 
parameter space than for the case of counting fractional 
number of such events in separate channels. The above 
is due in part because every SUSY event that passes the 
trigger has P™ ISS , and so one maximizes the signal events 
as opposed to obtaining a fraction of them. Further, the 
SM P™ lss falls off rapidly beyond the peak value coupled 
with the fact the P™ ISS is larger in SUSY extending out 
to momenta where the SM cannot produce a large num- 
ber of events. 

(ii) (Pj? lss ) Analysis: A remarkable signature emerges 
distinguishing the stau coannihilation region and the HB 
region if one analyzes Asusy for each parameter point 
as a function of (P™ ISS ) which is the mean P™ lss calcu- 
lated by averaging the P™ ss over the entire model event 
record. The above phenomenon is shown in Fig. (|TJ) (right 
panel). Here one finds that (P^ lss ) has a very wide range 
from 300 GeV to a TeV or more for the stau coannihi- 
lation region, while (Pj? lss ) for the HB region lies in a 



much narrower band centered around 350 GeV - a phe- 
nomenon which originates for parameter points on the 
Chargino Wall. Thus the (Pj? lss ) ranges in Fig. (JTJ) (right 
panel) can be viewed as one of the smoking gun signa- 
tures which can discriminate between the two mecha- 
nisms using LHC data. 

Although a quantitative analysis of (P™ ISS ) is rather 
complicated since it involves many particles and depends 
in part on post trigger level cuts, one can give a qual- 
itative picture of the disparity between the Pp lss on 
the Stau-Co and the HB regions by analyzing the decay 
chains of sparticles into their final products culminating 
into an odd number of LSPs (per sparticle decay chain) 
and the SM particles. Here one finds that often the spar- 
ticle decays on the Stau-Co involve two body decays. For 
the HB case, however, one finds that the sparticles pro- 
duced in pp collisions have typically a longer decay chain 
which depletes the P™ ISS in this case. 

We illustrate these features by analyzing the two spe- 
cific benchmarks given in the caption of Fig. Q (left 
panel). For the HB model point of Fig. Q (left panel) 
the following production cross sections are dominant: 
PP ^ (gg/xZxf/xtxf) at the level of (45,25,15)%. 
While squark production is highly suppressed (nig ~ 
622 GeV < m~ ~ r r G (1.2,1.7) TeV). One finds 
that the dominantly produced g decays via Br [g — > 
X°i + 9 + q] ~ 50% and Br[g -> xf + Q + ~ 50% 
with the LSP contributing only 10%. The reason for 
this largeness is because the on-shell decay of the gluino 
into qq is suppressed due to largeness of the squark 
masses (a phenomenon which typically holds for the 
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gluino decays on the HB). Further, the X2 an d xf P r0 ~ 
duced on HB have B_r[x£ -> X? + / + /] ~ 100% and 
Br[xf Xi + / + /'] ~ 100%. Thus the decay chain 
for sparticles produced on the HB tend to be longer re- 
sulting in reduced P™ ss Now, for the Stau-Co model 
point of Fig. fl]) (left panel) the leading SUSY production 
level cross sections are from (jjq, qq, gg) at the level of 
(41, 33, 7)%, with the corresponding 2 body decay modes 
Br[<7ft — > Xi + 9] ~ 100% (1st and 2nd generation), and 
Br[<7z, -> (xl-Xi ) + (<?,<?')] ~ (60,30)%. Since the decay 
chain for sparticles on the Stau-Co tend to be shorter 
the resulting P™ ISS is larger. In summary the main rea- 
son why the HB tends to give lower values of missing Pp 
relative the Stau-Co is simply due to the fact that on 
the HB, in order to get to the LSP from the dominant 
gluino production mechanism one usually needs at least 
2 successive 3 body decays, while on the Stau-Co the 
right-squarks (q"ii) from the first and second generations, 
which arc dominantly produced, each decays right into 
the LSP + quark. The above also holds more generally 
in that one finds that sparticles arising from the Stau-Co 
have much shorter decay chains resulting in fewer final 
particles and thus the missing energy can get large. Our 
more general results given here on a large spread in P™ lss 
agree with the first Ref. of [HI and with the analyses of 
the CMS and ATLAS collaborations 0, Ej^. Con- 
versely the models on the HB have longer decay chains 
with more final state particles and thus the missing en- 
ergy carried by the neutrals is depleted leading to missing 
P™ ss which is more SM like. This feature has also been 
discussed in the CMS analysis of [56| . 

(iii) Jet cuts, n* et : Another powerful signature for dis- 
criminating the Stau-Co and the HB regions is Asusy 
taken as a function of P|. miss and n* et , i.e., 

Asusy = Asusy (n jet > n* et , P™ ss ^ P} miss ), (2) 

where the * indicates a fixed cut value. Here features spe- 
cific to the coannihilation branch and to the HB emerge 
when p™ lss is fixed, and n* et is varied. This is shown 
in Fig.Q. In particular one finds that for the coanni- 
hilation branch there is an optimal n* et near 4 because 

the discovery limit criteria A|^ Y _Co /v / SM decreases as 
a function of increasing n* et and has a max near n jet ~ 4, 
while for the annihilations on the HB, specifically on the 
Wall, the situation is quite different, in that the larger 
the n*j et the larger is the value of A|^ SY /%/SM (where of 
course the SM is subject to the same * cuts). Thus, as 
the jet number n* et becomes large A™ SY sustains a much 

stronger signal than AfusV C °> and tllus A SUSY /\/SM is 
a strong discriminator between the Stau-Co and the HB 
regions. 

The above becomes very significant if the SUSY scale 
is high with the LSP mass lying in the several hundred 




FIG. 2: A discrimination of the two mechanisms for the sat- 
isfaction of the relic density in the early universe, i.e., stau 
coannihilation and annihilation on HB with assume d L HC 
luminosity of 10 fb _1 . Curves (connecting Nsvsy / VSM for 
discrete n* et ) correspond to various models with Mlsp < 
275 GeV for the Stau-Co and M LS p < 230 GeV for the 
HB. The ratio of A r ™sY tau-C °/v / SM is computed under the 
standard post-trigger level cuts (see also Fig.l) but with n* et 
taken as a variable. 

GeV range. Also, this type of large n-jet cut can deplete 
the leptonic signal, so a delicate balance of jet cuts 
is very important. In fact, if the SUSY signal is not 
highly leptonic, the analysis of the above type would 
be an efficient way to decipher new physics. Further, 
even if one has signatures with many leptons, the jet 
analysis will provide additional corroborating signatures 
for discovery and discrimination. 

(iv) Tagged b-jets: The utility of b-tagging for the HB 
region has previously been emphasized in [l5|, [l^, [52] • In 
Fig. ([3]) (upper panels) we give an analysis exhibiting how 
b-tagging provides a striking discrimination between the 
Stau-Co and HB regions where we plot A(nfe)/ v /SM(n&) 
as a function of the number of tagged b-jets (nb) and find 
this dependence to be drastically different for Stau-Co vs 
HB regions. In Fig. ©(lower panel) the fractional num- 
ber of events with 2b-jets vs the number of events with 
2b-jets is given, and again one sees a strong discrimina- 
tion between the parameter points in the Stau-Co region 
vs those in the HB region. In Fig. (J4j) we extend the anal- 
ysis to the 4b-jet mode and correlate this signature to 
events with two hadronically decaying tau jets and to 
events that do not contain tagged b-jets. 

As already discussed in (ii), on the HB gluino produc- 
tion is dominant in pp collisions at the LHC, and further 
one finds that the gluino decays dominantly into bb, i.e., 
Br[g -» x° + b + b] ~ 40% and Br[g -> xf + Kb) + i(t)\ - 
40%. Thus, the gluino 3 body decays are very rich in 
b quarks. Conversely for the Stau-Co model point of 
Fig.(l) (left panel) the pp production cross sections are 
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FIG. 3: (a): Top panels: N(nb) / ^SM(nb) vs nb for the Stau- 
Co and HB regions where N(nb) (SM(n&)) is the number of 
SUSY (SM) events that contain n b-tagged jets. A sharp 
discrimination between the Stau-Co and the HB by b-tagging 
is observed. The number n* et is fixed at 2. Here m-g < 1.1 
TeV. (b): Bottom panel: A plot of iV(26)/iVsuSY vs N(2b) 
where n^ et is fixed at 4. 



as follows: (gq) ~ 41%, (qq) ~ 33%, (gg) ~ 7% (as al- 
ready noted in (ii)). Further, the gluino has only a small 
branching ratio into bb in this case via g — > bb. Including 
the production cross sections for gg and the branching 
ratios, we find that overall the bb production on Stau-Co 
is smaller relative to that on the HB. More generally the 
analysis of Figs. [(JT])- (J2J)] shows that the LHC signatures 
arising from the Stau-Co are easily distinguishable from 
those arising from the HB. 

(v) N SVSY (leptons/jets) and (P™ iss ) vs af p : Signif- 
icant additional information regarding the coannihila- 
tion region and the HB region can be obtained by an 
analysis of the number of SUSY events at the LHC vs 
the spin-independent neutralino-proton cross section a^} p 
along with the current limits from the direct detection 
of dark matter. We give an illustration of the above 
in Fig. ([5]). The analysis of Fig. (right panel) shows 
that the coannihilation and the HB regions are well sep- 
arated in the space spanned by the trileptonic signature 
3L (L=e,fi) and <r^ p . Further one observes that in Fig. ([5|) 
(left panel) the parameter points in the HB region in the 
(P™ ss ) — cr^p plot are clustered together in a ball shaped 
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FIG. 4: N(2t) (the number of events with two hadroni- 
cally decaying r-jets) vs N(4b) (the number of events with 4 
tagged-b jets) (upper panel). A similar plot with N(2r) re- 
placed by N(0b) (the number of events with no tagged b-jets) 
(lower panel) . The VSM values in each channel are indicated 
by 'X' on the plots. Signatures arising from Stau-Co and from 
the HB regions are clearly discriminated. 



region and well separated from points in the Stau-Co re- 
gion which lie on a slope again providing a strong dis- 
crimination between the Stau-Co and the HB regions. 

In the above we have analyzed in detail how the LHC 
data can allow one to discover if the mechanism for the 
origin of dark matter in the early universe arises in the 
HB region or in the Stau-Co region. This opens up the is- 
sue if the above type analysis can be done more generally 
to identify the dominant mechanism for the generation of 
dark matter in the early universe using LHC data. Below 
we discuss briefly how one may extend the analysis to the 
stop coannihilation region and the A-pole region. Thus as 
already discussed in the introduction, in addition to the 
Stau-Co region there is also a stop coannihilation (Stop- 
Co) region where the relic density constraint is satisfied. 
However, regions of the parameter space which give rise 
to Stop-Co have signatures which are highly non-leptonic 
relative to those of the Stau-Co and of the HB [lj, 
Thus correlations such as OL + jets vs 1L+ jets allow 
one to distinguish Stop-Co regions from others as dis- 
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FIG. 5: Right panel: An exhibition of the trileptonic signal 
vs er^p. Points in the vertical region to the right constitute 
the Chargino Wall. Left panel: an exhibition of (P™ ISS ) vs 
o-^p. The cluster of points at the end to the right consti- 
tute the Chargino Ball. The CDMS/XelO constraints HH and 
constraints expected from SuperCDMS[K| are also shown. A 
clear discrimination of Stau-Co and HB can be seen in these 
plots. 

cussed in the analysis of [li], Such correlations if 

observed, would be a good indication of stop coannihi- 
lation as the origin of dark matter. We note in passing 
that the LSP is mostly a Bino in this case which sup- 
presses the scalar cross sections in the direct detection 
of dark matter [3l|. Next we discuss the A- funnel region. 
Here the relic density is satisfied because the LSP mass 
is nearly half the CP odd Higgs mass. The analysis of 
parameter points which cluster near the pole region do 
not have the same NLSP for all parameter points unlike 
the case of HB where xt is the NLSP, or the Stau-Co 
region where fi is the NLSP. However, in the pole region 
the NLSP could be Xi or ^i (*i is seldom seen in the 



pole region). Thus the A-po\e region can give mixed sig- 
natures, sometimes characteristic of HB and sometimes 
characteristic of the Stau-Co. To firmly establish the pole 
region one would need a global analysis with many sig- 
natures which would give a determination of the Higgs 
A mass and the mass of the LSP. A similar situation 
holds for other isolated regions of the parameter space 
which cannot be classified in the above categories which 
satisfy the relic density constraints. Here also one would 
need a global analysis on the signatures to identify the 
mechanism regarding the origin of dark matter. 

Conclusions: The analysis presented here shows that 
with sufficient LHC data one can discriminate between 
the Stau-Co and the HB regions regarding the origin of 
dark matter in the early universe in the mSUGRA model 
using lepton, jet and missing energy signatures. We dis- 
cussed several smoking gun signatures for such a discrim- 
ination. It was also shown that further discrimination is 
possible by combining LHC data with the limits on <r~£ 
from the direct detection of dark matter. An analysis in 
a similar spirit is given in [39| . Our analysis utilizes sim- 
ulations based on CMS detector specifications and spe- 
cific final state signatures in probing the origin of dark 
matter, and we do so over a large portion of the param- 
eter space. To our knowledge there are no analyses in 
the literature which have tried to utilize the expected 
LHC signatures to probe the origin of dark matter along 
the lines discussed here. While the analysis presented 
in our work illustrates our main points in the mSUGRA 
model, similar analyses along these lines should be pur- 
sued for other models of soft breaking including string 
and D brane models. Specifically, it would be interesting 
to analyze what the LHC can tell us about the origin 
of dark matter in cases where one has departures from 
mSUGRA including the case of a non-thermal relic such 
as a Wino dominated neutralino. 
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